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The bcc-based Fei-2,AL exhibit interesting magnetic and anomalous structural properties as a 
function of composition and sample processing conditions arising from thermal or off-stoichiometric 
chemical disorder, and, although well studied, these properties are not understood. In stoichiometric 
B2 FeAl, including the effects of partial long-range order (i.e., thermal antisites), we find the observed 
paramagnetic response (with non-zero local moments), in contrast to past investigations which find 
ferromagnetism based on local density approximation (LDA) to density functional theory or which 
find a non-magnetic state from LDA-I-U, both of which are inconsistent with experiment. Moreover, 
from this magneto-chemical coupling, we are able to determine the origins of the lattice constant 
anomalies found in Fei-^AIa, for x~ 25 — 50, as observed from various processing conditions. 

Although bcc-based iron-aluminum (Fei-ajAlj,) aUoys have been investigated extensively over the years, the complex 
(even anomalous) structural and magnetic properties as a function of composition and state of chemical order are 
not well understood. In fact, while most experimental investigations observe Curie- Weiss-like paramagnetic (PM) 
responses with small (~ 0.3^b) effective moments on FeiiSi^, electronic-structure calculations for perfectly ordered 
B2 (/3-CuZn) phase obtain a ferromagnetic (FM) state with ^ 0.7/is^*^*^i2i2i. In addition, experimentally in the 
iron-rich Fe-Al system, anomalies are observed in the average lattice constant that vary with processing conditions, 
i.e., whether the samples were quenched, annealed, and/or cold-workediS. The important fact is that the processing 
route used in preparing samples (stoichiometric or not) produces lattice defects that can have a significant effect on 
the mechanical, structural and magnetic properties of Fe-Al as well as on the intermetallics as a whole. Furthermore, 
there are competing structures: at 28% (38%) Al at 825 K (475 K), for example, the B2 phase undergoes a secondary- 
ordering transition to DO3 (FesAl or Heusler) phase. And, due to kinetic limitations, the phase diagram is known 
only for T > 475 K. To be able to compare to experiment and to gain an understanding of the anomalies found in this 
system will require a first principles method capable of including disorder and temperature effects on equal footing 
with the electronic structure. 

Here, using first-principles methods, we investigate the structural and magnetic properties of bcc-based Fe-Al, 
including the A2 (bcc disordered), B2, and DO3 phases, as a function of compositionally- and thermally-induced 
disorder (antisites). In B2 FeAl, we predict that the PM state competes with FM state when the state of partial long- 
range order is taken into account. In particular, at T= K, the PM and FM states are degenerate in energy at ~ 70% 
of perfect long-range order but are near-degenerate (< 0.25 mRy/atom) over the range of 50 — 90% of order (partial 
order compatible with experiment); whereas FM state is ~ 0.5 mRy/atom lower than the PM (and nonmagnetic 
(NM)) state for perfectly ordered B2, in agreement with past theoretical work. More generally, this magneto-chemical 
coupling leads to anomalies in the average lattice constant in ¥ei-x^\x for x ~ 0.25 — 0.5 at.%Al, reproducing what 
has been observed for annealed, quenched and cold-worked samples, but also never explained or found from theory. 
Our results provide a simple and physically reasonable understanding of the, heretofore, anomalous properties found 
in bcc-based Fe-Al. 



I. DISCUSSION 

Ordered intermetallics constitute an important class of high-temperature structural materials. One of the most 
common crystals structures in the intermetallics is the B2 structure, which because of its simple lattice structure makes 
it an ideal model for studying a variety of physical phenomena found in these systemsii. The physical properties of 
Fe-Al are quite sensitive to extrinsic defectsiSii^ and thermomechanical historjii2ii^ii^ii&. In fact, B2 FeAl is not fully 
ordered due to various lattice defects, such as vacancies and antisites, and they are thought to play a major role in the 
mechanical and magnetic behavior of this system. However, the precise concentration of the different lattice defects 
is not well known. Furthermore, unlike the strongly ordered intermetallics such as B2 NiAl, which forms triple defect 
structures consisting of two vacancies on the transition metal sublattice and an antisite defect on the Al sublattice, 
the possible defect structures for the less strongly ordered B2 FeAl are much more complicated, as the formation of 
both vacancies and antisite defects on the Fe site are thermodynamically stableii. 
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As mentioned, using the local density approximation (LDA) to density functional theory (DFT)i&, electronic- 
structure calculations on ideal B2-FeAl obtain a FM state with /i ^ 0.7 ng^i^t^tli^^ in contradiction to the Curie- 
Weiss-like PM response found in experiment. Calculations based on the so-called disordered local moment (DLM) 
paramagnetic state have also been suggested, but this approach by itself was unable to explain experiment®. Calcula- 
tions using generalized gradient approximation (GGA) to DFT, which are known to improve the magnetic ground-state 
description of elemental bcc FcfiS, have been employed but do not yield FM DOa-FesAl or NM B2-FeAl as a ground- 
state^li. On the other hand, recent investigations^^-^^ have obtained a NM ground-state for ideal B2 FeAl within 
so-called LDA-I-U procedure, which is an extension of LDA that introduces correlation corrections due to the orbital 
dependence of the Coulomb interaction. The LDA-I-U finds a NM ground-state for B2-FeAl for a particular range of 
the empirically chosen parameter U. However, once again, it is in contradiction to experiment, and does not confirm 
the LDA-I-U picture. Furthermore, none of these investigations address any other experimental observation, such as, 
the important composition-dependent structural anomalies strongly affected by processing. 

In alloys, both short-range order and partial long-range order (LRO), and associated magnetism, are strongly 
affected by the processing temperatures, whether its an anneal or quench, which freezes in chemical disorder at 
that temperature (e.g. < Tcurie < Tprocessing < Torder-disorder)- For example, DOa-FesAl crystal structure has 
^ 8% site disorder according to Bradley and Jay2^, which was as later confirmedS^^. Thus, the overall processing 
procedure affects greatly the materials properties and, therefore, the characterization data obtained from samples 
that are not completely ordered (due to kinetic limitations). The aforementioned DFT calculations provide results for 
perfectly ordered states or off-stoichiometric disorder corresponding to T=0 K, which may not be directly related to 
the experimentally assessed data. However, Johnson et. al.^^ have shown that including the thermally-induced partial 
LRO effects allows quantitative comparisons to characterization experiments, and explains discrepancy of T = K 
DFT results. Clearly an investigation of the effects of (thermally-induced or off-stoichiometry) partial order on the 
magnetic and structural properties in the B2 FeAl system is in order. 

II. COMPUTATIONAL DETAILS 

For our electronic-structure calculations, we use the Green's function based, multiple-scattering approach of Kor- 
ringa, Kohn and Rostoker (KKR)2L2£. We use the local density approximation to the exchange-correlation potential 
and energy as parameterized by Von-Barth-Hedin2S. In addition, all calculations are performed using the atomic 
sphere approximation (ASA) with equal size Fe and Al atomic spheres^. 

By using the KKR method, we can include chemical and magnetic disorder in the electronic structure and energetics 
via the coherent potential approximation''-'^ ''^ (CPA) modified to incorporate improved metallic screening due to 
charge-correlations arising from the local chemical environment''''. The KKR-CPA density- functional theory and 
total energy formalism provides reliable and well-documented energetics and structural-related parameters, from fully 
disordered to ordered configurations. 

We describe the PM state by the disorder local moment (DLM) approximation^^iSi^, where magnetic short-range 
order is ignored, yet the site-dependent average effect of magnetic orientational disorder is self-consistently included, 
whether chemically ordered (partially or fully) or disordered. Thus, our implementation of the KKR-CPA is capable 
of treating simultaneously chemical disorder, site defects and magnetic (orientational) disorder all on equal footing 
with the electronic structure (e.g., Slater-Pauling curved, or the INVAR''®-''^). 

For all calculations we use a four-atoms unit cell (see Fig. ^ based on an underlying bcc crystal lattice that 
inherently contains the DO3, B2 and A2 (bcc disordered) structures as a function of long-range order or composition 
because each sublattice is composed locally of Fei^^-^^x- Using Fig. ^it is easy to define order parameters needed. For 
example, referenced to the Al sublattice (i =A1), the B2 ordering at x=l/2 can be defined in terms of a temperature- 
dependent LRO parameter ?7(T): c^^; = i(l + 7j{T)) and cf^^i = i(l - r]{T)), giving the fully disordered A2 
(ordered B2) lattice for 77=0 (1). This describes a static compositional modulation of the A2 phase by a wavevector 
ko=27r(l,l,l)/a to yield the partially-ordered B2 FeAl ("pseudo-FeAl" ), see, e.g., Ref. E 

Characterization is performed on samples at finite temperatures with < '>]{T) < 1. The temperature dependence 
of ?y(T) can only be obtained from a thermodynamic calculation or measurement. In Fe-Al, ?7(T) varies continuously in 
the solid phase (for DO3 with arbitrary x, it can be done via two LRO parameters for the two operative wavevectors, 
i.e, (111) and (^^^)^''). Hence, a Landau expansion of the free energy difference in terms of F(7y) (relative to r]=0 
state) for a given fixed-sized unit cell, DOa-type in our case, is written as 

AF^ir^) = F(2)^-(0)?72+0(^4) 

where odd powers of rj are zero due to the translational symmetry of the A2 (ji — 0) phase, and AS is an entropy differ- 
ence for the magnetic state cr =(DLM, FM, or NM). We note that the largest contribution to the entropy is the point 
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Figure 1: The four-atom cell used to study cubic DO3. For basis atoms 1-^, the translational vectors in lattice units of 
Ti = (101), T2=(011), T3 = (001) generate the standard 16-atom (cubic) DO3 cell. In ideal DO3 Fe3Al, Fe (Al) occupies sites 
1-3 (site 4)', whereas in B2 FeAl, Fe (Al) atoms occupy odd (even) numbered site. Sites 1 and 3 are always equivalent. 

entropy and this cancels exactly within a common unit cell; also, when the nearest-neighbor (chemical) environments 
are approximately the same (as roughly is the case here), the nearest- neighbor pair entropy cancels exactly. Therefore, 
the temperature effects are primarily due to antisite disorder, and (for fixed composition) can be accounted for by 
analyzing total energy differences. (See Ref. |2^for how such calculations are related to characterization experiments.) 

In addition to calculations of B2 FeAl versus 77, fully ordered FeAl and FcsAl, we also consider Fei-ajA^ with Al 
concentration between 25% and 50%: fully-disordered and off-stoichiometric-disordered B2 and DO3 with and without 
chemical disorder on the Al sublattices (marked 2 and 4 in Fig- Specifically, the following cases of constituent 
concentrations are investigated: 



(i) For 0.25 < a; < 0.50: = = 0, cf ' = 1, = 4a; - 1; 

(ii) For 0.25 < a; < 0.50: cf = cf = 0, cf = cf = 2a;; 

(iii) For 0.25 < a; < 0.25-ha;o/2: cf = 1 - xq, cf = c^' = cf = (4a; - (1 -a;o))/3; 

(iv) For 0.25 + xo/2 < a; < 0.40: cf = l- xq, cf=cf=xo, cf={Ax - 1 - xq)- 

(v) For 0.25 <x< 0.25 + a;o/2: cf = cf = (4a; - (1 - a;o))/3, cf = cf = 0.5 - xq/S; 

(vi) For 0.25 + a;o/2 < a; < 0.50: cf = cf = xq, cf = cf = 2x - xq. 



Cases (i) and (ii) have no disorder on sublattices 1 and 3 and ideal (no LRO) DO3 and B2 compositions for off- 
stoichiometric Yei-^^^x, (iii)-(vi) are used to simulate possible compositions with LRO, i.e. thermodynamic chemical 
disorder. Concentrations of the constituents for sublattices 1 and 3 are assumed to be the same both in DO3 and 
B2, and the Al content is not more than xq (see Fig. QJcaption). The value of xo=0.09 is chosen somewhat arbitrarily 
for illustrative purpose, although it is comparable with the disorder reported experimentally from DOa-type Fe-Al in 
Refi^; for B2 FeAl it corresponds = 1 - 2xo = 0.82. 



As mentioned, experimentally B2 FeAl is found to be paramagnetic, with finite local moments. Yet, in agreement 
with previous theoretical work^i^iSiLS our calculations find that the ideal B2 FeAl is ferromagnetic. At the theoretical 
equilibrium lattice parameter a = 5.34 Bohr, we obtain a magnetic moment of w 0.72/xb for Fe and — 0.03/xb for Al, 
again consistent with the published data. However, such calculations ignore the fact that B2 FeAl is not fully ordered 
at any finite temperature. 

Therefore, in Fig. |21 we show the stability of FM, NM and DLM partially ordered phases of FeAl alloy versus 
long-range order, -q, relative to that of the fully-disordered FM state. 

The NM-FM energy difference decreases monotonically as 77 — *■ 1 (increasing order) from 4.8 mRy/atom in the A2 
state to 0.5 mRy/atom in the fully ordered B2 structure. For 77 = (A2) the DLM-FM energy difference, /\E^^^~^^ , 
is ~ 3.8 mRy/atom (a^*^ = 5.45 Bohr) with ^f^f = 1.83/xb and a^f^^^ = 1-32^5, respectively. We find that the 
FM state is nominally always the lowest energy state for any value of ry (^/^jt^dlm-fm ^ q 25 niRy/atom, i.e., less 
than thermal energies, for 0.9 > 7] > 0.5), except at 77 0.7 (1 — 7;^ = 0.5) where the FM state is degenerate with 
the DLM state. For this nearly ordered B2 structure, 77 — 0.7 corresponds to a 15% concentration of Al on the Fe 
rich sublattice, compatible with ~ 19% from experimen1j^ii£. Both the FM and DLM states have nearly the same 
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Figure 2: The DLM (circles), FM (squares) and NM (triangles) (a) total energies relative to the A2 FM (mRy/atom) and (b) 
Fe magnetic moments (hb) versus partial order 1 — 77^, with filled symbols the Al-rich sites (i.e., the Fe anti-sites). 

lattice parameter (5.36 Bohr). The local Fe anti-sites moments are finite for any partial long-range order, see Figure 
121 Whereas Fe atoms on Al-rich sites have nearly the same magnetic moments (~ 2.1/xb) for both FM and DLM 
states, the Fe moments for the other sites only exist in the FM state where they are relatively small. The average 
PM local moments are comparable to the observed Curie- Weiss moments ~ 0.3/iB at 77 = 0.7. 

It is clear that antisite defects play a crucial role in determining the magnetic configuration. However, in the DLM 
state at T] < 0.3 (see Fig. ^p), the Fe magnetic moments exist in both high- and low-spin states, similar to bcc Fe; 
while in the fully-ordered B2-phase (77 = 1), the DLM state is equivalent to the NM state, where all local moments 
are equal to zero (quenched). As an aside, we note that the PM state (stabilized, as identified here, by the antisite 
disorder) may benefit from further stabilization afforded by the LDA+U procedureSLSS, but it is not the origin for 
the PM effect, as is now clear. 

The key point is that the FM and DLM states can be degenerate or nearly degenerate in the partially ordered 
equiatomic FeAl state. Furthermore, Monte Carlo simulation^ of A2 FeAl (i.e. ry = 0) also obtain a paramagnetic 
state at room temperature. The important quantity in understanding non-zero temperature phenomena is the free 
energy and the disordered B2-DLM magnetic phase has a larger configurational entropy contribution to the free energy 
than the B2-FM ordered structure over a wide range of 77. Thus, partial long range chemical order is stable, making 
the paramagnetic DLM state the stable configuration at finite temperature, even with only a small degree of disorder. 

We note that the theoretical DLM lattice constant is about 3% below the experimental value (w 5.50 Bohr). 
However, such a deviation is comparable with pure bcc-Fe, as found in other calculations. The calculated DLM 
magnetic moment on an antisite Fe atom (i.e. on Al-rich sites) is close to the value of 2.2 /is observed by Parthasarathi 
and Beck^ for samples quenched in cold water. In addition, Bogner et al^ have extracted a moment of ~ 3.43/iB from 
Langevin-type behavior of magnetization-field dependence. Bogner et al. hypothesized that such a moment could 
result from a cluster formed by nine nearest Fe atoms, which would correspond to a partially-ordered alloy with a 
concentration of Fe atoms on Al sites of ~ 18%, close to the 19.12% found in experiment. Again, it is similar to the 
disorder we find needed for degeneracy in our calculations, i.e., ~ 15% with r] = \ — 2(0.15) 0.7. 
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Figure 3: Reduced lattice parameter of Fei-xA\x versus x. Theory results are for off-stoichiometric DO3 (triangles) and B2 
(circles), or the A2 phases (stars). Theoretical ao (5.320 Bohr) is the lattice constant of the A2 at x=0.25. Filled symbols 
(triangles/circles) are for DO3/B2 with additional thermal (antisite) disorder, see text. Experimental data is the dashed (dash- 
dotted) line for as-deformed samples (quenched from 1273 K)^", and ao is that observed lattice constant at x=0.25 for deformed 
sample. Quenching from 523 K yields similar results to high-T quench data. 



IV. STRUCTURAL ANOMALIES IN FE-RICH FE-AL 



Lattice constant anomalies have been experimentally observed in the Fe-rich region of the Fe-Al phase diagram. The 
lattice constant dependence on concentration is strongly dependent on the quenching or annealing procedurei^i used 
in synthesizing the material (dashed line in Fig. i.e. thermal antisites and vacancies. For low Al concentration the 
lattice constant has a nearly linear dependence on concentration. However, between 25%-40% the lattice parameter 
exhibits non-monotonic behavior with a maximum lattice constant at a; < 0.3 (for deformed samples at a; 0.34) and 
a minimum between 0.34 < a; < 0.43 with increasing Al content. Of course, the exact behavior strongly depends on 
the type of quenching procedure used during processingiS. On the other hand, the experimental data do not provide 
quantitative information on (chemical and magnetic) sample-dependent disorder in the observeci^^ partially ordered 
DO3 and B2 structures, which have chemical disorder on each independent sublattice (estimated to about 8% near 
X = 0.25) and partial long-range order. It is assumed that the formation of the "pseudo-FeAl" type of ordering is 
associated with the structural anomaliesiS. 

First, to address the anomalies, we consider the "ideal" DO3: case (i) - with disorder only on the Al sublattice 4', 
and ideal B2: case (ii) - with disorder on the Al sublattices 2 and 4- 

Our DO3 calculations exhibit linear growth of the lattice constant with increasing concentration up to x 0.29, 
while the average magnetic moment (per atom) /i decreases from 1.42/^^ at x = 0.25 to 1.23/iB. At x ~ 0.29 two 
degenerate solutions exist for the FM state. The low-spin solution has both a smaller /i ^ 0.54/LtB and a smaller 
lattice constant (see Fig. O. However, the largest discrepancy is for the Fe magnetic moments on sites 1 and 3 
where the low-spin solution has a significantly smaller magnetic moment on these sites (~ 0.2/ib) as compared to 
high-spin solution (~ 1.52/xb) on the same sites. Above x=0.29, the low-spin solution is energetically more favorable; 
at elevated temperatures a smooth transition from one solution to the second solution should be expected due to the 
possible coexistence of both solutions in a range of Al concentration. 

In addition we note that the calculated lattice parameter of the fully-disordered A2 phase is almost linear for 
concentrations of Al up to 34% which agrees well with experimental observations for deformed samples quenched at 
high temperatures, see Fig. O For the A2 phase the energy is about 7 mRy larger than "ideal" DO3 at x — 0.25, see 
Fig. 21 The respective energy difference begins to grow rapidly with increasing Al concentration at x 0.3: from 8.4 
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Figure 4: Formation energy of ferromagnetic A2, DO3 and B2 Fei-^AL versus concentration of Al; symbols are chosen as in 
Fig. El 

mRy aXx = 0.29 to 22.6 mRy at x = 0.5. 

Consideration of the "ideal" B2-disordered ¥ei-x^x reveals that the lattice parameter is practically indistinguish- 
able from the low-spin solution in DO3 calculations, however, with magnetic moments that are slightly higher than 
the ones found in DO3. The difference in the total energy between pseudo-FePA and "ideal" DO3 declines rapidly, see 
Fig. 21 from 2.9 niRy/atom at x=0.25 to 0.5 mRy/atom at x=0.40 and finally to zero at x=0.5, where the chemical 
structures are equivalent. Taking into account that pseudo-¥eK\ configurational entropy contribution to free energy 
is larger than in partially-ordered DO3 for the same composition, it is reasonable to expect that at finite temperature 
off-stoichiometric DO3 and B2 phases could coexist for some values of x. For high enough Al content DO3 would 
eventually be replaced by the B2 phase as the ground state where again the actual value for the Al concentration 
would depend on both the temperature and processing procedure. 

To illustrate the importance of the thermodynamical chemical disorder we performed calculations for both DO3 
(iii,iv) and B2 (v,vi) phases. Figure O shows a stronger lattice constant dependence on concentration for B2 than for 
the "high-spin" DO3 solution (i.e. at a; < 0.3). For larger Al content the influence on B2 and DO3 is comparable. In 
addition, the linear dependence of the formation energy on the Al concentration is virtually unaffected by the chemical 
disorder (excluding the uniform shift), see Fig. 01 for x > 0.34. However, for smaller Al concentrations this is not the 
case. This is due to the dependence of the sublattice composition on Al content in (iii)-(vi). The additional effect of 
partial-order on top of off-stoichiometric disorder leads roughly to a plateau in the formation energy, especially for 
DO3 (see Fig. ^J, which is compatible with the shift off stoichiomctry for the maximum in the B2-DO3 transition 
temperature. 

Generally, the magneto-chemical coupling effect is a ubiquitous phenomena observed in many materials. It has been 
used also to explain quantitatively, e.g., ordering characterization data™ and INVAR lattice anomaly in Fe-Ni^. 



V. SUMMARY 



The structural and magnetic anomalies observed in bcc-based Fe-Al have remained unexplained experimentally and 
theoretically for over 50 years. As characterization data and their interpretation depend upon the sample preparation 
and processing, the analysis of their properties based on simulation require the inclusion of these important thermal 
effects. We presented a theoretical study that incorporates magnetic, off-stoichiometric and thermodynamic chemical 
disorder on equal footing with the electronic structure, exemplifying a quantitative means for studying materials 
processed at high temperatures where point defects (such as anti-sites) play an important role in determining the 
properties. We find that partial long-range order due to thermal (i.e. antisite) disorder yields paramagnetic B2 
FeAl (with finite moments) that is stable with respect to the ferromagnetic state, in contrast to past theory but 
in agreement with experiment. The partial order required also agrees with experiment. In addition, we showed 
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that this same magneto-chemical couphng produces structural anomalies observed versus Al content. Our results 
provide a clear and simple explanation behind the observed paramagnetism and structural anomalies in Fe-Al. We 
expect that this type of long-range order effect is responsible for the observed FesAl lattice constant dependence on 
annealing temperature'*^, as well as the paramagnetism (or spin-glass behavior*'*) in B2-Fei_a;Al2;, both of which may 
be investigated using the approach employed here for Feo.sAlo.s. 
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